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Abstract 

In the past few years the B-factories became unexpected players in the scalar 
mesons business: in order to access the CP violation effects, it is necessary to 
handle the dynamics of the strong interaction between the final state hadrons. 
A number of large statistics studies heavy flavor decays involving a scalar com- 
ponent have been performed recently by Belle and BaBar, who have joined 
CLEOc, BES, E791 and FOCUS in the effort to understand the physics of the 
scalar mesons. In this talk, the most recent results from these experiments will 
be reviewed, with emphasis on the low energy Kn system and on the f0(1370). 



1 Introduction 



The identification of the scalar mesons is a long standing problem. There are 
too many candidates with mass bellow 2 GeV/c^, although some states, like 



the K and the /o(1370), are still controversial. In addition to the regular qq 
mesons, the observed spectrum may contain other types of states, like glueballs, 
hybrids, tetraquarks or molecules. 

There are difficulties from both experimental and theoretical points of 
view. Scalars are, in general, broad overlapping states. Since there is no spin, 
they decay isotropically. In scattering experiments the production rates for low 
mass states, like the a and k, are suppressed by the Adler zeroes. At higher 
masses, disentangling the broad, spinless states, like the /o(1370), from the 
smooth background is complicated by the interference with other scalars and 
with the continuum, which distorts the line shapes. 

The precise determination of pole positions and couplings to specific 
modes of all existing scalar particles is an essential step towards the identi- 
fication of the genuine qq states. Such an ambitious task could not be accom- 
plished by one single type of experiment. One has to look at the scalars problem 
from all possible perspectives, exploring the different constraints imposed by 
different reactions. 

In the past six or seven years there has been lots of new results on scalar 
mesons from heavy fiavor (HF) decays to light quarks (LQ), exploring the 
unique features of these processes. In this paper we will focus on two crucial 
problems: the nature of the k meson and the existence of the /o(1370). More 
specifically, we will present the latest results from r lepton and three-body 
decays of D and B mesons to light hadrons. 

The paper is organized as follows. In Section 1 we discuss why HF decays 
are a very useful tool for the study of scalar mesons. We will also discuss some 
aspects of the basic analysis techniques in HF decays. In Section 2 we discuss 
the situation of the low mass Kw spectrum. In Section 3 we discuss the issue 
of the /o(1370), in the light of hadronic D and B decays. The last Section 
contains a summary and conclusions. 

2 Heavy flavor to light quarks 

There are unique features that make decays of D and B mesons to light hadrons 
very suitable for the study of scalar mesons. These particles are copiously 
produced in D and B decays, especially when there is a pair of identical particles 
in the final state. The B-factories already liavc^ very large c;harm samples, with 
a high degree of purity. Soon there will be also large samples oi B ^ hih2hs 



decays {hi = tt, K) from the LHC expcrinicmts. With these high quahty data, 
the TTTT and Ktt spectra can be accessed continuously, starting from threshold, 
and covering the entire elastic range. 

Another appealing feature, especially in D decays, is the close connection 
between the quark content of the initial state and the observed resonances. In 
the decay of a D meson the weak decay of the c quark is embedded in a strongly 
interacting system that leads to the final state liadrons. However, if one goes 
through the PDG listing, one realizes that, in spite of the complexity of the 
D decay, nearly the entire hadronic and semileptonic rates can be described 
by a rather simple scheme, in terms of trcc-lcvel valence quark diagrams, and 
the regular qq mesons from the Constituent Quark Model. The dominant 
amplitudes in D decays are the external (spectator) and internal W radiation. 
The 'final state valence quarks' result from transition c s{d) + ud{s) plus 
the 'spectator' q. These quarks define the possible intermediate; qq states. This 
simple picture works very well for intermediate states having cither a vector, 
an axial-vector or a tensor resonance. If one excludes the scalar mesons, in 
D decays there is nothing else than the members of the usual qq nonets of 
the Constituent Quark Model. In other words, D decays can be seen as a 'gg 
filter': if a resonance is observed in D decays, then it is very likely to be a 
qq meson. One can expect that this holds also for the scalar resonances, so D 
decays would also provide clues about the nature of these mesons. 

Semileptonic decays of the type D — > h\h'2lv, and hadronic decays of the 
T lepton, r — > hih'2V, provide a particularly clean environment for the study 
of the scalar mesons, since there is no strong interaction between the 
pair and the leptons. However, the /11/12 system is dominantly in P-wave in 
both cases. The S-wave contribution is typically less than 10%, so very large 
samples are required. An additional difficulty is the fact that the neutrinos are 
not reconstructed, so the background level in these decays is relatively large. 

Hadronic decays of D mesons, on the other hand, are much easier to be 
reconstructed. In some final states the S-wave component is largely dominant, 
like in the D+ — > ii'~7r+7r+ and in the D+ — > 7r~7r+7r+ decays. Background 
levels arc typically of the order of a few percent. The problem here is how to 
disentangle the desired information. The final states are strongly interacting 
three-body systems, with a complex and unknown production dynamics. The 
pure /ii/i2 is certainly the main ingredient, but there is no direct route to extract 



it in a model independent way. Approximations, and, therefore, interpretation 
of the results, are unavoidable, unfortunately. 

Most of the existing data on HF LQ come from hadronic three-body 
D decays. The event distribution in the Dalitz plot is given by, 

^ ■M{Sa,St)\^ (1) 
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where Al is the mass of the decaying particle and Sa, st are the two-body 
invariant masses squared. The phase space density is constant, so any structure 
in the Dalitz plot reflects directly the dynamics of the decay. 

The analysis technique of such decays is by now standard ^ . The de- 
cay amplitude is written as a coherent sum of phcnomcnological amplitudes 
corresponding to the possible intermediate states, 
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= Ecfee^''=A^,^f = fk^ BWk (3) 

In the above equations fj^ and fj^ are form factors, with L being the 
orbital angular momentum at the D or at the resonance decay vertex; is 
a function accounting for the angular momentum conservation, and BWk is a 
relativistic Breit-Wigner function (with an energy dependent width, in general). 
The complex coefficients CkC^^'' are usually the fit parameters. 

Analyzes differ by the way the S-wave, A^, is modeled. There are three 
basic approaches. 

The most used model is the so called isobar model 12 © , in which the 
S-wave is represented by a sum, of Breit-Wigner functions and a nonresonant 
amplitude, with unknown complex coefficients. The nonresonant amplitude is 
assumed to be uniform in D decays, but varies across the Dalitz plot in the 
case of B decays. This model is simple and intuitive, but there are well known 
conceptual problems with this approach . As we will see, it provides an 
effective description of the data but, in some cases, the physical interpretation 
of the results is rather ambiguous. 



Another approach is the K-matrix model . This is a very sophisticate 
tool, but it is based on a very strong assumption: in the D decay, the reso- 
nant hih2 system and the bachelor particle recoil against each other without 
any interaction. This would greatly simplify the problem, although one must 
acknowledge the lack of experimental evidence supporting this assumption. If 
the three-body rescattering is negligible, the evolution of the /ii/i2 system must 
be the same as in /11/12 elastic scattering. In the K-matrix approach the S-wave 
amplitude is, therefore, fixed by other type of reactions, whereas the produc- 
tion of the hih2 pair is an unknown function to be determined by the fit. This 
universality of the S-wave is often seen as the most appealing feature of the 
K-matrix approach. We should not forget that the same constraint should 
be applied to all partial waves, which, unfortunately, is not the case of the 
analyzes published so far. In general, good fits are achieved only if an extra 
energy-dependent phase is added. In the framework of the K-matrix model, 
the origin of such a phase is attributed to the production of the /ii/i2 pair. 
This is simply a matter of interpretation, since this phase could also be due to 
the rescattering of the final state particles. There is no way, with the existing 
measurements, to distinguish between these two effects. 

Finally, there is the P WA approach'*^^ . Here no assumption is made about 
the content of the S-wave, which is treated as a generic complex function of 
the /11/12 mass, .4° = a{s)e'^'^^^\ The /11/12 spectrum is divided in bins. At 
each bin edge the amplitude is defined by two real constants, and 4>k, which 
are fit parameters. The amplitude at any value of the /11/12 mass is given by 
a polynomial interpolation. This method relies on a precise modeling of the 
P- and D-waves. The problem resides, once more, in the interpretation of the 
results. The measured S-wave phase (f>{s) includes any rescattering/production 
effects, which should be disentangled in order to determine the 'bare' hih2 
amplitude. 

One last remark is in order. A consequence of representing the decay 
amplitude by a coherent sum of amplitudes is that the decay fractions do not 
add to 100%. Different amplitudes populate the same region of the phase 
space, so they are expected to interfere. The amplitudes have phases that vary 
across the Dalitz plot. The interference can be destructive in some regions of 
the phase space and constructive in other ones. One should be very careful, 
though, when the sum of fractions largely exceeds 100%. In almost all cases this 



large interference occur between the broad amplitudes in the S-wave, involving, 
in general, the nonresonant amplitude. This is a symptom of a poor modeling 
of the S-wave. In other words, one may find a mathematical solution to the 
fit problem, but with a misleading physical interpretation. In the following 
Sections we will see several instances of this problem. 

3 The low energy Kn spectrum: the k problem 

There are two crucial questions to be answered in the low mass Kit spectrum. 
The neutral n has been observed by different experiments in several types of 
heavy flavor decays^ El El)^ go far evidence for its charged partner has 
been scarce and controversial. This casts doubts about the k being a regular 
/ = 1/2 qq state. The other issue is the pole position of the n. There is no 
data on Kit Kit elastic scattering bellow 825 GeV/c^. Determination of the 
K pole position ^ must rely on extrapolations of the existing data. New data 
that fill the existing gap would be highly desirable. This section is devoted to 
these two issues. 

3.1 The D° K+K-tt" decay - BaBar 

In a recent analysis of the decay 1)° — > K~^K~7t'^, performed by the BaBar 
collaboration ® , the issue of the charged k has been addressed. The data 
sample has IIK signal events with 98% purity. The Dalitz plot of this decay 
is shown in Fig. [1] Resonances can be formed in all three axes (the third axis, 
K~K"^ , is along the diagonal, starting at high K^tt^ mass) of the Dalitz plot. 
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Figure 1: Dalitz plot of the — > K^tt^t:^ decay ^ . The narrow band at the 
upper edge of the Dalitz plot correspond to the ^tt" channel. 



The main diagrams for this decay are shown in Fig. [2l There is a well 
defined pattern in D decays that proceed via the external ly-radiation am- 
plitude (Fig. [2]-a): due to the V-A nature of the weak interactions, the W 
couples much more strongly to a vector or an axial-vector meson than to a 
pseudoscalar. According to this pattern, one expects a dominant contribution 
from the K* (892)^ K~ channel, compared to K* (892)^ K'^ . Important contri- 
butions from the K* {892)~ K'^ and (pir''^ modes are also expected. We can see 
in Fig. [T] clear structures corresponding to these resonances. 

The Dalitz plot of Fig. [T] was fitted using three different models for the 
K^tt'^ S-wave. In the first fit, the S-wave was represented by the usual isobar 
model - a constant nonresonant amplitude plus two Breit-Wigner functions for 
the K='=7r° and i^o (1430)='=7r° modes - with parameters taken from E79l'31). The 
second fit used the E791 PWA S-wave® measured from the 13+ — ^ _Rr~7r+7r+ 
decay. Finally, the third fit used the LASS 1=1/2 K^n^ S-wave ampli- 
tude^I^. The isobar model has the smaller fit probability [x^ prob. < 5%). 
The E791 PWA S-wave provides a good description of the data (x^ prob. = 



23%), but the best fit was obtained using the LASS 7=1/2 S-wave. Results 
of the fits with LASS 7=1/2 S-wave are summarized in Table [T] Note that 
in Model II the exclusion of the K* {14:1Q)K amplitudes - a 5% contribution 
in Model I - has a minor impact on the other P-wave components, but causes 
a dramatic change in the S-wave fraction. Moreover, in Model II the sum 
of the decay fractions largely exceeds 100%, indicating the existence of large 
destructive interference effects. 



Figure 2: Valence quark diagrams leading to the K^K^n'^ final state: external 
(left) and internal (right) W radiation. The resonances that can be formed by 
the 'final state' quarks are indicated. 

The fit with the LASS amplitude is much better than the one having 
explicitly the amplitude. Note that this result does not rule out the charged 
K, since its pole can be found in LASS data It simply says that the isobar 
representation of the K^tt^ S-wave in the 7?+ K^tt^tt^ decay, from E79ll3 
and FOCUS '^15 analysis, is not a good model for the TC+tt^ S-wave amplitude 
in the L>" ^ K+K-7r° decay. In both E791 and FOCUS analyzes of the 
7?+ — ^ 7r~7r+7r+ decay, the LASS 7 = 1/2 S-wave amplitude fails to provide 
a good description of the data, whereas a very good fit was achieved with the 
isobar model. The origin of this discrepancy is yet to be understood. With 
a larger D'^ K^K^tt*^ sample, a model independent measurement of the 
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Table 1: Decay fractions (%) of the — > K^K vr'' decay. Results are from 
fits using the 1—1/2 LASS S-wave amplitude. 



mode 


model I 


model II 


K*{892)+K- 


45.2±0.9 


44.4±0.9 


K*iU10)+K- 


3.7±1.5 




K+n'^iS) 


16.3±0.1 


71.1±4.2 




19.3±0.7 


19.4±0.7 


/o (980)^0 


6.7±1.8 


10.5±1.4 


K*{892)-K+ 


16.0±0.9 


15.9±0.9 


K*{1A10)-K+ 


2. Til. 5 




prob. 


62% 


47% 



S-wave amplitude could be performed and directly compared to the results of 
the E791 MIPWA®. 

3.2 The K^TT^Vr decay - Belle 

In hadronic decays of the r lepton, r —> /ii/i2'^r, the hih2 system is not affected 
by strong interaction with the leptonic current. This would be as close as one 
could get to the ft.1/12 /11/12 elastic scattering using HF decays. The Belle 
Collaboration published recently t^S' a study of the decay r~ li^'K~Vr- A 
sample with 53K signal events was selected from the reaction e^e~ — > r+r^. 
The decays to a muon plus two neutrinos, so the signature of the event is 
one lepton recoiling against three charged prongs in the opposite hemisphere. 
With a total of three missing neutrinos, it is very difficult to reconstruct the 
event topology. Immediate consequences are a high background level (~20% 
in this analysis) and the lack of an angular analysis. 

The angular distribution of the helicity angle was used in FOCUS studyllS 
of the decay. The helicity angle is formed by the mo- 

mentum and the line of flight of the Z?+, measured in the K~tt'^ rest frame. 
In this decay there is a 5% scalar component. The line shape of the K~tt^ 
spectrum, which is dominated by the if* (892)°, is not sensitive to the different 
models for the S-wave component. However, the different possibilities - a com- 
plex constant, a Breit-Wigner function, or the LASS 7=1/2 S-wave amplitude 
- lead to different angular distributions, so one can explore this feature in order 
to understand the nature of the scalar component. 



In the case of the t~ ^ K n~i^r decay, the K*{892)~ alone is not enough 
to describe the K tt" spectrum, as we can see in Fig. [3]-a. There is an excess 
of events in the lower and in the upper part of the spectrum. The K tt~ 
spectrum was fit with two classes of models : the dominant K*{892)~ plus the 
charged k and one high mass state (either a vector or a scalar resonance); the 
i^*(892)- plus the LASS 7=1/2 K-Tr+ S-wave amplitude. The best fit was 
obtained with the K*{892)^ plus a pure S-wave, that is, a sum of the k and 
the ifo(1430)". The confidence level of this fit is 41% (Fig. Hb). The fit using 
the LASS amplitude, on the other hand, yielded a C.L. of only 10^^. 




Figure 3: a) The K n mass spectrum'^^ (points with error bars). The his- 
togram represents the A'*(892)^ contribution, b) The k'^'k~ mass spectrum 
with the different background components and the result of the best fit superim- 
posed. 



The result of Belle analysis is in conflict with those of the BaBar study 
of the K^K^TT^ decay. One could expect that, in the absence of strong 

interactions, the K tt" S-wave phase would match that of LASS, whereas in 
the case of D° K^K~tt^ the rescattering would cause deviations from the 
elastic phase. This is a rather intriguing result. In both cases the statistics is 
limited to a few thousands of events, though. The angular analysis on a larger 
^ K Tx^ Vr sample may confirm the resonant behavior of the S-wave at low 



mass. In this case we would have a compelling evidence of the charged k(800). 
3.3 The D+ -> /f-7r+7r+ decay - FOCUS 

The — » K^TT^TT^ is a golden mode for the study of scalar mesons. The S- 
wave accounts for over 60% of the decay rate. In addition, the large branching 
fraction (9.5%), combined with the long £)+ lifetime and a final state having 
only charged tracks (the odd charge track being always a kaon), makes the 
Z?+ K~TT^TT^ the easiest charm decay mode to be reconstructed. Using this 
decay one can measure the K~'!t~^ S-wave amplitude near threshold, filling the 
existing gap in LASS data. 

FOCUS has published recently a study of this decay using the K-matrix 
approach 'liJ' . A good fit was obtained combining the 1—1/2 and 1 = 3/2 
LASS S-wave phases with an additional energy dependent phase. This extra 
phase was interpreted as being originated from the production dynamics of the 
Ktt system, but it could also be attributed to the rescattering of the final state 
particles. One important aspect should be stressed. The 1 — 3/2 component 
is purely nonresonant, whereas all resonances are in the 7=1/2 amplitude. 
The latter has also a nonresonant background. The fractions of each isospin 
component in the K-matrix fit are (207±28)% and (40± 10)% for the 7=1/2 
and 7 = 3/2, respectively. The total S-wave fraction amounts to (83±2)%. Here 
is another instance of large interference between broad components within the 
S-wave. 

The PWA approach for the 7C^7r+ S-wave was used recently by CLEOcttl 
and FOCUS l^'. In this new FOCUS study, the D+ R- 71+11+ Dalitz plot 
is fitted using also the isobar model for the K~'k+ S-wave. This is a work in 
progress, based on a sample of 93K signal events and with 98% purity. From 
here to the end of this Section we will discuss the preliminary results this 
FOCUS PWA/isobar Dalitz plot analysis. 

In the isobar fit, the S-wave model was that of E791 a sum of an 
uniform nonresonant amplitude, plus two Breit-Wigner functions for the K7r+ 
and 7^0(1430)71+ modes. In order to make a direct comparison with E791 
and CLEOc, the S-wave Breit-Wigner functions are multiplied by the same 
Gaussian form factors, fo — ^'"d/i^ Masses and widths of the scalar 
resonances are fit parameters. 

The fit fractions and resonance parameters are shown in Table [51 Results 



Table 2: Decay fractions (%) of the — > K^tt^'k^ decay. Results are from 
fits using the isobar model for the S-wave amplitude. The masses and widths 
of scalar resonances are in units of MeV/c^ . 



mode 


E791 


CLEOc 


FOCUS(a) 


FOCUS(b) 


FOCUS(c) 


7?'(892)07r+ 


12.3±1.4 


11.2±1.4 


11.3±0.3 


11.7±0.3 


11.2±0.3 


i?*(1410)"7r+ 
7?*(1680)°7r+ 






1.2±0.3 


1.1±0.3 


1.3±0.3 


2.5±0.8 


1.4±0.2 


3.3±0.3 


2.7±0.3 


3.8±0.3 


7?2(1430)7r+ 
;^o(1430)7r+ 


0.5±0.2 


0.4±0.4 


0.20±0.05 


0.20±0.05 


0.20±0.05 


12.5±1.4 


10.5±1.3 


16.8±0.8 


14.3±0.7 


18.7±1.2 


K(800)7r+ 


47.8±13.2 


31.2±3.6 


43.5±4.5 


71.3±5.5 


22.3±3.2 


nonresonant 


10.4±1.4 


13.0±7.3 


14.3±3.0 


7.5±3.1 


31.6±4.5 


k(800) mass 


797±48 


805±11 


837±12 


829±14 


867±14 


k(800) width 


410±97 


453±21 


443±21 


433±18 


485±27 


l?o(1430) mass 


1461±3 


1459±5 


1466±4 


1468±4 


1466±4 


l?o(1430) width 


169±5 


175±16 


193±7 


193±7 


192±7 



from the three experiments arc in reasonable agreement, except for the S-wave 
fr'action in CLEOc analysis, which is a bit smaller than E791 and FOCUS. But 
the most remarkable feature in Tabled is shown in the last three columns. The 
isobar fit to FOCUS data was repeated with a little variation of the S-wave 
parameterization. In FOCUS(b) the value of the scalar form factor parameter, 
ro, was set to 6 GeV-^ (instead of 5 GcV^i in FOCUS(a)). In FOCUS(c) the 
value of was set to zero, which is equivalent of having no scalar form factor. 
All the three FOCUS fits have an equally good confidence level. We observe a 
dramatic change in the S-wave composition, while the P- and D- waves remain 
unaltered. The nonresonant fraction varies by a factor of four. 

This instability can be readily explained by the interference between the 
broad k Breit-Wigner and the uniform nonresonant component. In Fig. |3] 
the effect of the scalar form factor is illustrated. Fig. |3]-a shows the modulus 
squared of the k Breit-Wigner. The Ktt mass dependence of the width, in 
the denominator of the Breit-Wigner, shifts the maximum of the function to a 
value bellow the Kit threshold. The introduction of the Gaussian form factor 
modifies the line shape near threshold (Fig. [4]-b). The resulting k amplitude 
has a more reasonable behavior. The same effect can be achieved if, instead of 
multiplying the k Breit-Wigner by the Gaussian form factor, we add the correct 



amount of a constant nonresonant amplitude, with the right phase difference. 
This is shown in Fig. O-c. 

The case of the Z?+ K^tt^tt^ is didactic: the isobar model yields a 
good fit, but it fails to provide an unambiguous physical picture of the decay 
dynamics. The model for the S-wave is clearly inadequate. We need to go 
beyond the isobar model in order to understand the composition of the S-wave. 

The £>+ i<'~7r+7r+ Dalitz plot is also fitted using the PWA method. 
The Ktt spectrum is divided in forty equally spaced intervals. At the edge of 
each interval the S-wave amplitude is defined by two fit parameters, Ao{s — 
Sk) = OfcC*"^''. The value of the S-wave amplitude at any point in the Ktt 
spectrum is given by a spline interpolation of these forty points. The S-wave is 
determined by an iterative procedure. A first fit is performed fixing the P-wave 
to that of the isobar analysis. The S-wave is determined (80 fit parameters). 
A second fit is performed fixing now the S-wave and varying the P-wave. The 
process is repeated until it converges. 




Kit mass (GeV/c') Kti mass (GeV/c^) Kn mass (GeV/c') 

Figure 4: Left: The modulus squared of the k Breit-Wigner, without the scalar 
form factor; center: the modulus squared of the k Breit- Wigner multiplied by 
the scalar form factor; right: the modulus squared of the sum of the k Breit- 
Wigner, without the scalar form factor, and a complex constant. 

Fig. [5] shows the FOCUS PWA S-wave phase, (j){s), as a function of 
the Ktt mass squared. A 80° overall phase was added to 4>{s) for a better 
comparison with the 1=1/2 and / = 3/2 phases from LASS. There is a clear 
discrepancy between the S-wave phase measured by FOCUS and that of LASS 
1=1/2 amplitude. The agreement between FOCUS and LASS cannot be 



achieved even combining both isospin phases. An extra phase, which depends 
smoothly on the Kn mass, is necessary for the matching of the two amphtudes. 
The origin of such a phase is unclear. 
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Figure 5: The phase of the S-wave Kir amplitude, from the 
decay, as a function of the Kit mass squared. The phase was shifted up by 80 
degrees, in order to make a comparison with the LASS 1—1/2 S-wave phase, 
Si/2- The LASS L — 3/2 S-wave phase, (S3/2, is also shown. 



The PWA is as close as one can get to a model independent method, 
since no hypothesis is assumed for the S-wave. But the method has some 
technical limitations arising from the large number of free parameters. One 
limitation is the existence of multiple solutions. With so many parameters, the 
fit has freedom to cure eventual problems with the P-wave model. A precise 
representation of the P-wave is mandatory, otherwise some 'leakage' into the 
S-wave would be unavoidable. 

The interpretation of the results is not trivial. The amplitude measured 



with the PWA method includes not only the Ktt dynamics, but also any pos- 
sible contribution from the production of the Ktt pair and rescattering of the 
final state particles. Some urgent input from theory is required in order to 
disentangle these effects. 

4 Is there an /o(1370)? 

In the low energy tttt spectrum the a pole is now well established, but the 
situation between 1-2 GeV/c^ is still controversial. There may be three neigh- 
bor scalar states, namely the /o(1370), the /o(1500) and the /o(1710). The 
/o(1500) is a narrow, well established state, with mass, width and couplings 
know to a good degree of accuracy 113' , but the uncertainty on the /o(1370) 
parameters is very large: 1200 < toq < 1500 GeV/c^ and 200 < Tq < 500 
GeV/c^. The existence of this state is often questioned. 

The /o(1370) has been observed mostly in central production and pp 
annihilation (see, for instance D. Bugg's recent review on this state SZ])) BES 
has also reported on this state from the decay J/^p ipirir This state is 

very difficult to be detected because it is broad and very close to the /o(1500). 
Its line shape may also be sensitive to the opening of the a Api channel. 

An interesting and related issue: the mass of the lightest scalar glueball is 
expected to be around 1.5 Gev/c'^. It is widely accepted that the three scalars 
would mixed states, having both qq and gg components in their wave functions. 

Heavy flavor decays are particularly useful here, since in these decays the 
qq component of the intermediate resonances are probed. 

4.1 The D+ 7r-7r+7r+ decay - FOCUS and E791 

The 7r^7r+7r+ decay is particularly suited to the study of scalar mesons. 

In this decay the S-wave component amounts to over 80% of the decay rate. 
The dominant diagram is shown in Fig. [6]-a, with a small contribution from 
the annihilation diagram (Fig. [6]-b). The /o(980)7r+ should be the dominant 
mode. The /o(980) has a large coupling to K^K^ and, therefore, it must have 
a strong ss component in its wave function. An important contribution from 
the channel /2(1270)7r+ is also expected. 

FOCUS has collected a sample of 1400 ^ 7r^7r+7r+ signal events. 
The Dalitz plot, shown in Fig. [3 was fitted with the K- matrix approach It 



was also fitted using the isobar model, but the results of this analysis were not 
published. 

There are two remarkable features in Fig. [71 a very clear structure at 1 
(GeV/c^)^, corresponding to the /o(980); the concentration of events near the 
border, at m^^ - 2 (GeV/c^)^, which is due to the /2(1270), p{U50)° and to 
a scalar state, which we refer to as fo{X). The Dalitz plot was fit using 
the isobar model for the S-wave, which has three components: the /o(980)7r+, 
the /o(X)7r+ and the nonresonant modes. The mass and width of the fo{X) 
are fit parameters. 




Figure 6: The dominant diagrams leading to the Df ir tt'^'tt^ decay. 

The fit fractions, as well as the mass and width of the fo{X), are shown in 
Table [3] (FOCUS errors are statistical only). FOCUS has about twice as many 
events as E791. In both analysis the /o(980)7r+ is the dominant component, 
followed by the /o(^)7r+. The BES CoUaboration^ISI) 

measured the mass and 

width of the /o(1370) in the J/tp ^ (pirir decay, obtaining toq = (1350 ± 50) 
GeV/c^ and Tq = (265 ± 40) GeV/c^. Comparing these values to the ones of 
the fo{X), we conclude that the state observed in the Z?+ 7r^7r+7r+ is not 
the same as the one observed by BES. 



Several variations of the S-wave model were tested. The crTr"*" mode was 
added to the S-wave, but its contribution is consistent with zero. A fit in- 
cluding the /o(f370)7r+ (with BES parameters) was performed, yielding a null 
contribution of this mode. In both FOCUS and E791, only one fo{X) is neces- 
sary to describe the data, and this state is consistent with being the /o(1500). 
The decay fraction of the fo{X) is very large. Given the diagram of Fig. [5]-a, 
one may conclude that this state has a significant ss component in its wave 
function. In this case we may also expect a large fration of this mode in the 
Df 7r^7r+Tr+ decay. 
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Figure 7: The Dalitz plot of the Df n 7r+7r+ decay, from FOCUS. 



4.2 The D+ K-K+tt+ decay - BaBar 

The Z?+ 7r^7r+7r+ and the Df K^K^n'^ decays share the same dom- 
inant diagram (Fig. [6]-a). The latter mode can also proceed via the internal 
PF-radiation amplitude (Fig. [5]-b). A dominant contribution from the <?!>7r"'", 
/o(980)7r+ and /o(X)7r+ modes is expected, but there should also be a large 
'K*{^92fK+ component. 



Table 3: Decay fractions (%) of the Df — > tt 7r+7r+ decay. Results are from 
fits with the isobar model for the tt^tt^ S-wave amplitude. 



mode 


FOCUS 


E79ll2Q|) 


/o(980)7r+ 


76.9±4.9 


56.5±5.9 


fo{X)7T+ 


23.3±0.5 


32.4±7.9 


nonresonant 


13.2±5.7 


1±2 


p(770)O7r+ 


1.2±0.1 


5.8±4.4 


p(1450)°7r+ 


4.0±1.0 


4.4±2.1 


/2(1270)^+ 


9.7±1.4 


19.7±3.4 


mo(/o(X)) (GeV/c^) 
ro(/o(X)) (GeV/c2) 


1.476±5.7 
119±18 


1.434±18 
173±32 



BaBar collected a very large (lOOK signal events) and clean (95% purity) 
sample ^21]) of the decay Df K^n^ . The Dalitz plot is shown in Fig. [9l 

BaBar is currently analyzing this data using the PWA method for the S-wave. 
Here we present preliminary results of the Dalitz plot analysis using the isobar 
model. 




Figure 8: The dominant diagrams for the decay — > K K~^t:^ . 

The fit result is shown in Table [H In Fig. [9] we have the Dalitz plot 
projections with the fit result superimposed. In the low mass region 

there could also be a contribution from ao(980)7r+, in addition to the /o(980)7r+ 



Table 4: Decay fractions (%) of the K K^tt^ decay, from a fit using 
the isobar model. 



mode 


fraction(%) 


/o(980)7r+ 
/o(1370)i^+ 
/o(1710)i^+ 
0(lO2O)7r+ 


35±14 
6.3±4.8 
2.0±1.0 
37.9±1.9 
48.7±1.6 



and the ^tt"''. In practice, is nearly impossible to fit the data with a model 
having, at the same time, these three amplitudes. The interference between 
them is very large, the coefficients become highly correlated and the individual 
fractions become too unstable. A stable fit is obtained with a model having 
only one of the two scalar amplitudes. The values reported here are from a fit 
with the /o(980)7r+. Note that there is still a large uncertainty in the /o(980)7r+ 
fraction. 

The most surprising result is the absence of the /o(X)7r+. We can see in 
Fig. [S]-b that there are very few events in the faiX) region. There is, on the 
other hand, a small excess of events next to the 0//o(98O) region, which is not 
well described by an uniform nonrcsonant amplitude. Instead, a scalar state 
was introduced. The fitted mass and width of this state are mp = (1.313 ± 
10 ± 114) GeV/c^ and Tq = (0.395 ± 8 ± 133) GeV/c^. The large errors reflect 
the sensitivity to the details of the S-wave parameterization. One could not 
really interpret this result as an indication of the /o(1370), since this is a very 
complicated region of the K^K~ spectrum. 

The absence of a fo{X) contribution in K^K^ti^ may indicate 

that the S-wave model used in the study of the _D+ 7r^7r+7r+ is not the most 
correct. BaBar is currently analyzing a sample of the decay D'^ 7r^7r+7r+ 
which is a factor of 10 larger than that of FOCUS. The most important result 
would be a simultaneous PWA measurement of the S-wave in both final states. 
The statistics is not a problem for the -D+ K^n^ , but, unfortunately, it 

is still a bit limited in the case of — » 7r^7r+7r+, even considering the BaBar 
sample. 
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Figure 9: a) The Dalitz plot of the Df K-K+tt+ decay ^ . Plots h) to 
d) show the projections of the Dalitz plot into the three axes (points with error 
bars) with the fit result superimposed (solid histograms). 



4.3 The D" k\+tt+ decay - Belle 

Belle and BaBar have collected very large samples of the K^tt+tt" 
decay. The Dalitz plot analysis performed by both experiments used the isobar 
model, and the results are in very good agreement. Here we will discuss the 
Belle analysis, based on a sample of 534K events with 98% purity. 

The diagrams for this decay are shown in Fig. 1101 The dominant ampli- 
tude should be the K*{892)~tt^ channel, with important contributions from 
the p{770)K mode and form the tt~tt~^ S-wave. The Dalitz plot, shown in Fig. 



1111 is very complex, since there could be resonances in all three axis. More- 
over, there is a small contribution from the doubly Cabibbo suppressed decay 
D° K'-'tt^tt^. In Fig. [TT]the label mt refer to the K 7t~ mass squared, if 
the parent is a D", or to the the K'^tt^ combination, in case the parent is a 




Figure 10: Diagrams for the K tt+tt decay. 



The fit model has 19 amplitudes leading to the K tt+tt" final state. The 
TT+TT^ S-wave contains four amplitudes: uK , the /o(980)ii' , fo{X)K and an 
extra scalar state, the . The fo{X) parameters were taken from E791 (see 
Table [3]), whereas the <72 parameters were determined by the fit. This extra 
CT2 amplitude was introduced to account for a structure near rrij^T^ ~1 GeV/c^, 
but it does not correspond to a real state. The parameters obtained by the fit 
are mo = (1.059 ± 6) GeV/c^ and Fq = 0.059 ± 10 GeV/c^. 

The dominant contribution is, as expected, the K*{892)^tt^ (62%), fol- 
lowed by the p°K° (21%) and the TT+TT S-wave (^^15%). The contribution 
of the fo{X) is small but significant (1.6%). No errors on the fractions were 
quoted. 
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Figure 11: The Dalitz plot of the decay ^ K tt+tt" . The label to^ 

refer to the "K^tt" mass squared, when the parent is a D'^ , and to the K'^tt^ 

— 

mass squared, when the parent is a D . 

In spite of the large number of amplitudes, including the extra '(T2', a 
good fit was not obtained. It is very hard to obtain a good C.L. in fits to very 
large samples. The goodness-of-fit is accessed by x^-like tests, in which the 
phase space is divided in bins of variable area, so that the number of events 
is similar in all bins. The isobar approach is, perhaps, too simplistic. With 
such a large sample, relatively small and localized deviations from the observed 
Dalitz plot distribution have, in general, large impact in the fit C.L. 

The goal of this analysis is to study the mixing phenomenon. The tech- 
nique is a time dependent Dalitz plot analysis. For this purpose, an effective 
representation of the data would suffice, according to the authors. Unfortu- 
nately the treatment given to the tt~^tt~ S-wave do not allow us to draw any 
conclusion. We cannot interpret the fraction attributed to the foiX) as an 
evidence for this state, since the tttt S-wave is not well understood. No attempt 
to measure its parameters was reported. A study of the D'^ K tt^tt" decay 



focused at the tttt S-wave and using the PWA method is in order. 



4.4 Charmless hadronic three-body decays of B mesons. 

Charmless hadronic three-body decays of B mesons are a very promising tool, 
but there is still a long way to go. The data samples resemble those of D 
mesons from the late 80's. There are two main problems: statistics is still 
limited and the background is still high. The nonresonant component is another 
problem. It is likely to be larger in B than in D decays. A constant nonresonant 
amplitude is the usual parameterization in the case of D decays, which may 
be a good approximation given the limited phase space. In B decays, however, 
the understanding of the nonresonant amplitude is a crucial issue , as one 
can already conclude from the existing data. 

There has been intense activity in this area, with many studies from 
the B-factories. Here we will focus on two analyzes from Belle, the i?+ 
i^T+TT+TT" and B° i^°7r+7r~ t^SI decays, and on two analysis by BaBar, 
the B+ K+K+K- ^ and B° ^ K°K+K- ^ decays. 

The B~^ K^TT^TT~ signal and the Dalitz plot are shown in Fig. [T^l 
Fig. [T21-b illustrates how large is the phase space and how much the action is 
concentrated near the border. 

In these decays the dominant mechanism for the b —^ s transition are 
penguin diagrams, shown in Fig. 1131 The diagram of Fig. 1131 a leads to final 
states having three kaons. The Kmr final states proceed via the diagram in Fig. 
[Hb. We expect dominant contributions from the ii'*(892)7r+, K^{U30)t:+, 
p{770)K and foK modes, in addition to the nonresonant component. We also 
expect the decay fractions in both i?° — * K'^tt^tt" and i?+ — > K^tt^tt" to be 
similar, since the replacement of the d by the u quark in Fig. 1131 b turns the 
B° to the B+ decay. 

All studies of charmless hadronic three-body decays of B mesons are per- 
formed with the isobar model. The nonresonant is parameterized by empirical 
formulae, with independent coefficients for each axis. In the two B Kirn 
analyzes by Belle the expression for the nonresonant amplitude is Anr = 
iie"'V(si) + a2e"'^/(s2)- In both studies the data is better described by a 
model having one tt+tt^ scalar state at m^jr ^ 1.5 GeV/c^. The results of 
Dalitz plot fit are in Tables [5] and [B) 
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Figure 12: a) The B+ K+tt+tt- signal from Belle\M . b) The B+ 
K^tt^t:^ Dalitz plot. 




Figure 13: a) Dominant diagram for B KKK decays b) Dominant diagram 
for B — > Kim decays. 

The decay fractions are similar in both decays, and correspond to the 
modes one expect from the diagram of Fig. [T51b. In both B^ and _B+ decays 



Table 5: Decay fractions (%) of the X+tt+tt decay, from a fit using 

the isobar model for the S-wave amplitude. 



mode 


fraction(%) 


K*{892)-Tr+ 


13.0±1.0 




65.5±4.5 


p{nofK+ 


7.9±1.0 


fo{9S0)K+ 


17.7±3.6 


MX)K+ 


4.1±0.9 


nonresonant 


34.0±2.7 



Table 6: Decay fractions (%) of the — > if "tt+tt decay, from a fit using the 
isobar model for the S-wave amplitude. 



mode 


fraction(%) 


K*{S92)-'K+ 


11.8±1.7 




64.8±7.8 


p{n{)fK+ 


12.9±2.0 


/o(980)if+ 


16.0±4.2 


.UX)K+ 


3.7±2.4 


nonresonant 


41.9±5.5 



there is a large interference between the ifo(1430)7r+ and the A'tt nonresonant 
component, which causes the fraction of the -ftTg (1430)7r'^ to be very high. The 
Kt: S-wave seems to be not well understood. The interference between the tttt 
S-wave and the corresponding nonresonant amplitude is small, though. The 
/o(^) state was represented by a Breit-Wigner function, whose parameters 
were determined by the data: toq = 1.449 ± 0.013 GeV/c^ and Tq = 0.126 ± 
0.025 GeV/c^. These values are in good agreement with the ones obtained 
from the D+ n-T:+n+ decay by FOCUS and E791. 

It is interesting to these results to those from B KKK . From the 
diagram in Fig. 1131 -a one expect significant contributions from the (pK and 
/o(980)if . The nonresonant amplitude is parameterized by empirical formulae 
similar to those used in Belle analysis. In the case of the B^ the nonresonant 
amplitude has three independent terms. Like in the B — > if tttt decays, a scalar 
K^K^ resonance was introduced, with mass and width determined by the fit. 

The K^K^ projections of the Dalitz plot are shown in Fig. [TH We 
see a clear bump next to the ((> peak, at ruKK ^^1.5 GeV/c^. The fit results 



are shown in Tables [7] and [S] The contribution of the (t)K decay is similar in 
both cases. As in the case of the Z?+ K^K^tt^ decay, the fraction of the 
fo{98Q)K mode suffers from large uncertainties. There is a small contribution 
from the fo{171Q)K , but only in the _B+ decay. The most striking features, 
though, are the very different K^K^ S-wave composition and the enormous 
interference between the K^K^ nonresonant term and what BaBar calls the 
Xo(1550) state. In the fit the sum of decay fractions amounts to 300%!. 

The Breit-Wigner parameters of the Xo(1550) state were determined in 
the i?+ analysis. In the analysis the Xo(1550) parameters were fixed at 
the values obtained in the analysis. These are: mo = 1.539 ± 0.020 and 
To = 0.257 ± 0.033. GeV/c^. The Xq{1550) parameters are very diflterent 
from the ones of the fo{X) state in i? — ^ Ktttt. However, before drawing any 
definitive conclusions, the K^K" S-wave is understood, as well as the role and 
form of the nonresonant amplitude. 




Table 7: Decay fractions (%) of the -B+ — » K K^K decay, from a fit using 
the isobar model for the S-wave amplitude. 



mode 


fraction(%) 


(t>K+ 


11.8±1.2 


/o(980)if+ 


19±8 


Xo{l550)K+ 


121±20 


/o(1710)X+ 


4.8±2.9 


nonresonant 


141±17 



Table 8: Decay fractions (%) of the — »• K^K^K decay, from a fit using 
the isobar model for the S-wave amplitude. 



mode 


fraction(%) 


/o(980)i^o 
Xo(1550)if° 
/o(1710)/^o 
nonresonant 


12.5±1.3 

40.2±9.6 
4.1±1.3 

112±15 



5 Summary and conclusions 

Two of the most challenging problems in the scalar mesons physics have been 
discussed from the point of view of heavy flavor decays. 

In the low energy Kir spectrum, the neutral k is now established. An 
analysis of the elastic scattering data revealed the position of the neutral n 
pole, in spite of the lack of data bellow 825 GeV/c^, and of the suppression 
of the amplitude due to the Adler zero. The existence of the n charged part- 
ners remain unsettled, though. These issues can be addressed by heavy flavor 
decays. The Ki^ S-wave amplitude was measured in the — > K~'!t^iT'^ de- 
cay, including the region bellow 825 GeV/c^ where LASS data starts. There is, 
however, one unavoidable task: to extract the K-k elastic scattering phase from 
the measured amplitude one has to handle other strong interaction effects. 

Evidences for a charged k are still scarce. The recent results from BaBar 
[D^ — *• K~K^'K^) and Belle (t~ K 'JT~i'r) are rather intriguing. In the r 
decay the Kir system is isolated from any other strong interaction. We would 
expect that the Ktt phase from this decay to match that of LASS, whereas 
in the case of the D° decay the three-body FSI, or a complex production 



amplitude, would cause some deviations from the pure elastic scattering phase. 
The experimental results, however, show exactly the opposite picture. More 
data and refined analysis techniques are clearly necessary. In the case of ther 
decay, the angular analysis is the crucial and missing piece. If such analysis 
is performed and confirm the resonant behavior at low Kn mass, the we will 
have a compelling evidence for the charged k. 

In the TTTT system, the existence of the /o(1370) has been also addressed 
in studies of HF decays. The analysis of the Df — > 7r~7r+7r+ and B — > Kirn 
decays show that only one scalar state with mass near 1.5 GeV/c^ is necessary 
to describe the data. The measured parameters of this foiX) state are consis- 
tent with those of the /o(1500). If we exclude the scalar mesons, in three body 
D decays the only intermediate states observed are those having a regular qq 
resonance. We could say that it is very likely that the /o(1500) is a genuine qq 
meson, or has, at least, a strong qq component in its wave function. On the 
other hand, the evidence for the /o(1370) in HF decays is weak and inconclu- 
sive. The puzzling fact, though, is that the fraction of the /o(1500)7r+ mode 
in the Df 7r~7r+7r+ decay is very large, suggests that the /o(1500) has a 
fairly large ss component and, therefore, a significant coupling to KK. But 
this is not true for the /o(1500), and no indication of this state was foimd in 
K~K~^n^ decay. One possible interpretation is that the S-wave model 
used in the Df -jt~w~^'jt~^ analysis is incomplete. There are also indications 
of a scalar state in B ^ KKK with mass near 1.5 GeV/c^, but no conclusions 
can be drawn before the nonresonant component is understood. 

We have seen that hadronic decays of heavy fiavor are a very rich envi- 
ronment for the study of the scalar mesons. Thanks to their unique features, 
the information provided by HF decays are complementary to the traditional 
hadronic collisions. One must also keep in mind that HF decays have been 
the only new data available in the past ten years, and this will be so until the 
commissioning of the new facilities. The B-factories already have very large 
and clean samples of Z) — > /ii/i2/i3 decays. As for the B — » hih2hz, data with 
equivalent quality will be available in a few years, from LHCb and the other 
LHC experiments. 

The existence of good data, however, is not enough. The experimentalists 
need to develop better analysis techniques, going beyond the isobar model. The 
limitations of the latter appear either when one moves to really high statistics. 



or when complex final states arc analyzed. Even with a better understanding 
of the S-wavc, it seems that a simple coherent sum of amplitudes may be 
too simplistic. The improvement of the analysis techniques does not depend 
only on the experimentalists creativity, but also on a deeper understanding of 
the decay dynamics, of the role of final state interactions, of the nonresonant 
amplitude, form factors and line shapes. This requires the urgent intervention 
of the theoreticians. 

In high energy physics most of the attentions are turned to the searches 
for new physics. There is a widespread belief that we are on the verge of 
great discoveries, that new particles are right at the corner. One of the most 
promising fields is the phenomenon of CP violation. The correct measurement 
of the CP violation effects, however, depends on the accurate understanding of 
the low energy strong interaction dynamics. That's where the flavor physics and 
the hadron physics communities meet. Even with somewhat limited analysis 
tools, there are plenty of good data from the B-factories that should be analyzed 
in a systematic way, with the focus on the physics of the light quarks and, in 
particular, of the scalar mesons. This is a very rich and challenging program 
that needs to be implemented. 
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